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ABSTRACT 


A flight simulation system with graphics and software on Silicon Graphics computer 
workstations has been installed in the Flight Vehicle Design Laboratory at Tuskegee 
University. The system has F-15E flight simulation software from NASA Dryden which uses 
the graphics of SGI flight simulation demos. On the system, thus installed, a study of pilot 
induced oscillations is planned for future work. Preliminary research is conducted by 
obtaining two sets of straight level flights with pilot in the loop. In one set of flights no 
additional delay is used between the stick input and the appearance of airplane response on 
the computer monitor. In another set of flights, a 500 ms additional delay is used. The 
flight data is analyzed to find cross correlations between deflections of control surfaces and 
response of the airplane. The pilot dynamics features depicted from cross correlations of 
straight level flights are discussed in this report. The correlations presented here will serve as 
reference material for the corresponding correlations in a future study of pitch attitude 
tracking tasks involving pilot induced oscillations. 
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INTRODUCTION 


McRuer and Graham (1981) have presented a survey of eighty years of flight control, and Ashkenaz 
(1984) has presented a survey of twenty five years of handling qualities research. From both survey 
presentations, it appears that prior to 1984, the cross correlations between pilot inputs and airplane 
response on a real time simulator had not been studied, but transfer functions representing human 
pilots had been proposed (see, for example, Washizu and Miyajima, 1965, and Smith, 1965), and the 
pilot induced oscillation (PIO) theory and resulting prediction technique had been developed (see 
Smith, 1977). The prediction techniques were tested on eighteen aircraft/flight control system landing 
configurations by Bjorkman (1986). Bjorkman indicated the need of more data and simulator studies 
to gain physical insights into PIO mechanization. From more recent works on PIO including the 
paper by Hess and Kalteis (1990), the need of more data is further desired. Cardullo et al, and 
Middendorf et al are some of the recent studies on the effect of simulator transport delay on pilot’s 
performance. In the present work a low cost flight simulator system with graphics and software on 
computer workstations has been installed. On the simulation system, a future study of PIO is 
proposed. As prepatory work for a future PIO study, the flight data has been obtained for two sets of 
flights with pilot in the loop. For conducting their flights, the pilots utilized the heads-up display on 
the computer monitor and operated a joystick to maintain a given altitude and a given heading. One 
set of flights is obtained without any additional delay between the stick input and the appearance of 
airplane response on the computer monitor. In another set of flights, a 500 ms additional delay is 
used. The flight data has been processed to find cross correlations between deflections of control 
surfaces and response of the airplane. It is anticipated that these cross correlations for straight level 
flights would serve as reference material to compare the respective cross correlations for future 
simulator studies of PIO. A portion of the work of this report was presented at an AIAA Southeastern 
Student Conference By Williams and Johnson (1997), and their paper is included as Appendix B in 
this report. While the test flights for the present study were being conducted, the raw data on some of 
them were displayed on the internet at http://silicon.tusk.edu/~jpwill. 


THE SIMULATION SYSTEM 

The simulation system is comprised of two SGI workstations and a BG System Joystick called a 
FlyBox. The two workstations are SGI IRIS with 2 CPUs and SGI Indigo 2 with one CPU; they are 
connected through ethemet. The simulation software, which represents the dynamics of flight of an 
F-15E jet fighter has been provided by NASA Dryden Flight Research Center. For the graphics of the 
system, the graphics of the SGI flight simulation demos has been integrated with the software of the 
airplane dynamics. Norlin (1995) has provided a description of the design of software, and its 
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simulation capabilities. According to Norlin, the F-15E real-time simulation is coded with a 
FORTRAN 77 shell and C support routines and it operates on a UNIX-based multiprocessor 
computer. The FORTRAN code includes the aircraft models, equations of motion, integration table 
look-ups, initialization and display generation routines. The C support routines are used for the 
graphical user interface, memory mapping, priority boosting, and interrupt handlers. In addition, C is 
used for the graphics and distributed system functions. The integration scheme is optimized for real- 
time operation. In this method, the derivatives are only calculated once for each frame, a weighted 
average of the midpoint and previous frame derivative is used to predict the derivative at the end of 
the frame. The equations of motion assume flat earth and six degrees of freedom. The airplane is 
assumed symmetrical about its own x-z plane. The simulation model calculations and integration of 
the equations are performed in the real-time loop. For pilot in the loop or hardware in the loop, the 
main real-time task is interrupt driven at the highest allowable system priority. The use of 3 CPUs 
ensures that timing constraints are adequately met. For their simulators, NASA Dryden prefer to use 
the name engineering simulators instead of production training simulators. The engineering 
simulators provide high fidelity simulation that is responsive to the needs of the researchers. The F- 
15E simulator used in the present work is an engineering simulator. It is called the F-15 ACTIVE 
which implies that it incorporates Advanced Control Technology for Integrated Vehicles. The code 
has the flexibility to program a desired transport delay between the stick input and the airplane 
response on the monitor. The use of UNIX “makefiles” to compile the simulation allows minor 
updates within about 5 minutes. The joystick has three different movements to deflect the three 
control surfaces: forward-rearward movement for elevator deflection, right-left movement for aileron, 
and clockwise-counterclockwise rotation about its own axis to simulate pedal movement for rudder 
deflection. It may be noted that the experienced pilots would prefer the use of rudder pedals instead 
of the existing stick movements to simulate the pedal movements, therefore, the absence of pedals 
may cause a lack of coordination between the rudder and ailerons in their test flights during 
simulation. In a simulated flight, the digital data on forty two different flight parameters can be 
recorded at every 25 ms time intervals. The flight parameters include altitude (H), time rate of change 
of altitude (HDOT), and heading (PSI) as well as stick deflection for elevator (DEP), stick deflection 
for aileron (DAP), and pedal movement for rudder (DRP). A pull on the stick or a positive value of 
DEP deflects the elevator upward which raises the nose of the airplane and increases its altitude (see, 
for example, Roskam 1995, p. 241). In the present report the DEP, DAP, and DRP data are essentially 
used to find cross correlations. It is understood that the deflections of elevator, aileron, and rudder are 
proportional to the values of DEP, DAP, and DRP, respectively. Therefore, for computing the 
dimensionless cross correlation coefficients, the values of DEP, DAP, and DRP are treated as relative 
values of deflections of elevator, ailerons, and rudder although the dimensional values of DEP, DAP, 
and DRP are not the dimensional values of deflections of the respective control surfaces. 
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FLIGHT DATA PROCESSING 


In the present study, all flights conducted on the simulator were straight level flights of three minute 
duration. For any flight conducted, the simulator can provide digital data for 42 flight parameters at 
every 25 ms time interval. For the present work, however, the digital data was obtained for only 
altitude (ft.), heading (deg.), time derivative of altitude (ft./sec.), and the deflections of elevator, 
rudder, and aileron. For the altitude, H represents its instantaneous value and h represents the 
fluctuating part of its instantaneous value. Therefore, 

h(t) = H(t) - avg. H 

where the avg. H is the average altitude of the three minute flight and t represents time. 

For the heading, X F represents its instantaneous value and \|/ represents the fluctuating part of its 
instantaneous value. Therefore, 

\|r(t) = 'P(t) - avg.'F 

Similar relations are used between instantaneous values, fluctuating values and average values of the 
control surface deflections. A is used for instantaneous value of control surface deflection, and 8 for 
the fluctuating part of the instantaneous value. For distinction between the three control surface 
deflections, the subscripts a and r represent aileron and rudder, respectively. No subscript is used on A 
or 8 to indicate elevator. The rms values of the above signals are denoted by, 8 \ h\ and v|/‘. The 
coefficient of cross correlations between 8 and h for an arbitrary time gap % is denoted as CORR (t) 
and it is defined as: 


CORR(t) = 


1 

T 8’ h 1 


t=T 

J 8(t) h(t + T) dt 

t=0 


The simulator provides data at every 25 ms interval, thus a 3 minute flight has 7201 points for each 
signal. The coefficient of cross correlation is therefore calculated as follows: 
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CORR(j) = 


1 

(7200 -j) 5' h’ 



To cover a time gap, x of 20 seconds, CORR(j) is calculated for j = 1, 2, ..., 801. For 100 second time 
gap j ranges from 1 to 4001. It is suggested that care be exercised in interpreting the long time gap 
correlations. For the correlation with zero time gap, the averaging is based on 7200 pts, with 20 
second time gap the averaging is based on 6400 points, and with 100 second time gap the averaging 
is based on 3200 points. The cross correlation coefficients between any other two signals are defined 
and calculated in the same manner. In the above formulas 5 and h may be replaced by any other two 
signals between which the cross correlation is desired. In the present work, the cross correlations are 
considered between 8 and h, 8 and dh/dt, 8 r and \}/, 8 a and \|/, 8 r and 8, 8 a and 8, 8 r and 8 a . 


THE PILOTS AND THE TEST FLIGHTS 

Four volunteer pilots conducted the required flights on the F-15E real-time simulation. Three of 
them were experienced licensed pilots and one had no prior experience of flying. In this report, the 
volunteer pilots are named PI, P2, P3, and P4. The ages, flying hours, and the license status of the 
pilots were PI: 69 yrs, 3200 hrs, SEL, MEL, INST; P2: 48 yrs, 100 hrs, SEL, MEL, INST; P3: 23 yrs, 
830 hrs, SEL, MEL, INST; and P4: 26 yrs, no flying experience. 

In every test flight, a pilot was required to maintain straight level flight for three minutes duration 
heading due North, with preadjusted throttle for flying at 450 knots indicated at 10,000 ft. altitude 
with autotrim OFF. During the flight, the pilot inputs were limited to the joystick movements only. 
The volunteer pilots were provided with the facility to train themselves on the simulator by 
conducting several three minute flights until they reached their own asymptotic performance in terms 
of (IT + 10\|f) a weighted sum of the rms values of fluctuations in altitude and heading in ft. and 
deg. respectively. Every pilot conducted two sets of straight level flights. One set of flights was 
obtained without any additional delay between the stick input and the appearance of airplane 
response on the computer monitor. In another set of flights, a 500 ms additional delay was used. 
Including the flights with the asymptotic performance, the pilot PI conducted 15 flights without 
delay and 15 flights with delay; P2 conducted 16 flights without delay and 21 flights with delay, P3 
conducted 22 flights without delay and 17 flights with delay; and P4 conducted 40 flights without 
delay and 25 flights with delay. From amongst the asymptotic performance flights, for each of the 
two kinds of flights, and for each pilot, based on the value of (h’ + lOx]/’), the best flight, the average 
flight, and the worst flight were picked up for further analysis. For all twenty four of these flights, 
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Williams and Johnson (1997) have reported average values and rms values of altitude and heading 
and 8-h correlation versus time gap graphs. 

A comparison of the average and rms values of altitude and heading for the four pilots indicated that 
the pilots P3 and P4 had remarkably small deviations, whereas pilots PI and P2 had rather large 
deviations from the desired values. It is believed that the cross correlations between any two randomly 
changing but related variables would be better representative of their mutual interdependence if their 
deviations from the desired quantities are rather small. Therefore the correlations for the flights of 
pilots PI and P2 are not reported here due to their large deviations in 8 and h. They may instead be 
seen in Williams and Johnson’s paper provided in Appendix B. In the main report, the discussions on 
correlations dwell upon five different flights and those are the best flights of pilots P3 and P4 without 
delay, the best flights of pilots P3 and P4 with delay, and the average flight of pilot P3 with delay. 
The average flight of pilot P3 with delay is included in the discussion because the best flight of P3 
showed an unusual correlation feature when compared with the other flights of pilots P3 and P4. The 
graphs of the raw data for these selected five flights are provided in Appendix A. 


FLIGHT TEST RESULTS 

For comparison of pilot’s performance in flights without transport delay and those with 500 ms 
delay, the rms values of altitude and elevator deflections for the best flights of the four pilots are 
shown in Table I. For both kinds of flights, pilots P3 and P4 have achieved better performance 
compared with pilots PI and P2. For the flights of pilots P3 and P4, the rms values of altitude and 
elevator deflection for flights with delay are approximately twice their respective values for flights 
without delay. It is obvious that the 500 ms transport delay has resulted in deteriorated pilot’s 
performance. With the existing related publications in view, this is not surprising. For example, 
Middendorf et al (1991) have noted that 300 ms transport delay resulted in significant deterioration 
of pilot’s performance. 

In a study of pilot dynamics and pilot’s choices of feedback alternatives for controlling altitude, the 
frequency response functions of pitch angle and altitude have been obtained and used by Goto and 
Matsuo (1988). In the present work, the cross correlations between the control surface deflections and 
airplane reponses are obtained to understand some aspects of the pilot dynamics. Let us start with the 
consideration of cross correlations between elevator deflection, 8 and airplane altitude, h. An upward 
instantaneous deflection of elevator from its average deflection is positive 8 and a downward 
deflection is negative 8, although the convention is to take downward elevator deflection as positive, 
for example Etkin and Reid (1996, p. 33). A higher instantaneous altitude compared with the average 
altitude of the three minute flight is positive h, and a lower than average is negative h. It is understood 
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that a positive change in 8 results in a positive change in h and a negative change in 8 results in a 
negative change in h. Allowing for a time lag between elevator deflection and altitude response, the 
two signals 8 and h should have a positive correlation. 

Figure 1 presents 8-h correlations for pilots P3 and P4 for a range of time gaps between zero and 20 
seconds. Figure 1 has the graphs for the five selected flights, their selection is indicated in the 
previous section. All the flights of Figure 1 have a common feature - that they do not have large 
deviations from the required altitude of 10,000 ft; the maximum deviation in average altitude is 6 ft. 
in the best flight of pilot P4. All the flights show an appreciable magnitude of negative 8-h 
correlation for zero time gap. This is expected because upon noticing positive h, the pilot moves the 
stick to deflect the elevator downward and vice versa. 

Let us consider the best flight without transport delay for pilot P3 as a typical one amongst the five 
flights of Figure 1 . In this flight the largest magnitude of negative correlation is reached when the 
time gap is 400 ms; to talk about it, let us call it the first time gap on the graph. The first time gap, 
perhaps, includes three stages of airplane altitude control: the pilot’s reaction time, the time gap 
between stick movement and the resulting elevator deflection, and the time gap between elevator 
deflection and affected altitude change. This consideration appears reasonable when it is seen with 
Roskam’s (1995, p.765) assumptions of 100 ms reaction time for test pilots and 120 to 200 ms 
reaction time for other pilots. The correlation stays negative up to 2700 ms. Let the second time gap 
on the correlation graph be the time gap between the largest negative magnitude of correlation and 
zero correlation, here it is 2700 - 400 = 2300 ms. The second time gap, perhaps, represents 
limitations on the rates of climb and descent; higher rates would lead to lower value of the second 
time gap. With transport delay, the best flight of pilot P3 has an exceptional feature that its’ first time 
gap is zero; this exceptional feature is not explored further. The raw data in Appendix A, however, 
shows no significant differences in the best and average flights of pilot P3 with delay; but the elevator 
deflections for the first 30 seconds in the flights with delay are different from the respective values in 
the flight without delay. Due to the exceptional nature of the best flight of pilot P3 with delay, the 
average flight of P3 is included in Figure 1 to be discussed as another typical flight with delay. Four 
of the five flights in Figure 1 may be treated as typical straight level flights; for them the first time 
gap ranges from 350 ms to 600 ms and the second time gap ranges from 2300 ms to 5700 ms. One 
of the reasons for large differences in the first and second time gap values for different flights may 
be that pilot’s observations and actions are continuous rather than discrete in nature; the pilot does 
not need to wait until elevator deflection or airplane response corresponds with certain stick 
movement. The largest magnitudes of correlation coefficients and the first and second time gaps are 
not significantly affected by the 500 ms transport delay. This is surprising, especially when we note 
that the rms values of fluctuations of altitude and elevator deflection in flights with transport delay are 
approximately twice the corresponding values in flights with no delay. The respective rms values are 
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provided in Figure 1. It is understood that a change in elevator deflection causes a pitching moment 
that brings about a rate of change of airplane altitude. Therefore, for a more clear view of certain 
features of the same flights, to substantiate the S-h correlations of Figure 1, the 8-hdot (hdot = dh/dt) 
correlations are provided in Figure 2. 

In Figure 2, all the flights, except the best flight with transport delay for pilot P3, show an 
appreciable negative 5-hdot correlation at zero time gap. This is expected because upon noticing 
positive hdot above the required 10,000 ft. altitude, the pilot moves the stick to deflect the elevator 
downward, and upon noticing negative hdot below 10,000 ft. he moves the stick to deflect the 
elevator upward. The negative values of 8-hdot correlations at zero time gap in Figure 2 are less in 
magnitude than the corresponding 8-h correlations in Figure 1. This is understood because upon 
noticing positive hdot for altitude below 10,000 ft. and negative hdot for altitude above 10,000 ft., 
the pilot does not need to change the elevator deflection. 

In the best flight of P3 without transport delay, picked up here as a typical flight, the correlation 
reaches zero when the time gap is 400 ms. For the same flight, this reinforces the observation that the 
largest magnitude of negative 8-h correlation occurred at the same time gap as seen in Figure 1 . In all 
the graphs of Figure 1 and Figure 2, it is noticeable that the 8-hdot correlation becomes zero at the 
same time gaps at which the 8-h correlations have local maximum or minimum values. 

In Figure 2, let us consider the two flights of pilot P4. For the 8-hdot correlations the number of zero 
crossings in 20 seconds for the flight with transport delay is twice that for the flight without transport 
delay. A similar feature is not found in the two kinds of flights for pilot P3. It may be reiterated here 
that pilot P3 is a certified and experienced pilot whereas the experience of pilot P4 is limited to flight 
simulators only. From Figure 1, we recall that the first time gap goes from zero to a value where the 
magnitude of negative correlation is largest and the second time gap goes from the end of the first 
gap to a value where the correlation is zero. Let the subsequent time gaps represent the subsequent 
zero crossings of the correlation starting from the first zero crossing. Perhaps the time gaps between 
the consecutive zero crossings would provide an indicator of the pilot’s performance, larger gaps 
representing better or more steady performance provided that the deviation in altitude stays small. 
Such an indication is jeopardized to a certain degree because stick signal in the flights under study is 
contaminated by electronic noise. For elevator and aileron deflections, the noise signal is estimated at 
3% of a typical stick movement. For rudder deflection, the signal is significantly larger than 3%. No 
attempt is made to suppress or filter out this noise because it requires the pilot to remain vigilant if he 
wishes to maintain the straight level flight. To find the subsequent zero crossings of the 8-h 
correlations, the graphs are obtained for time gaps from 0 to 100 seconds. For the same flights which 
are included in Figure 1, the 100 second graphs are shown in Figure 3. As suggested in the section on 
Flight Data Processing, care is needed in interpreting the long time gap correlations. In Figure 3, for 
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the time gap range from zero to 100 seconds, for pilot P3, a flight with transport delay has 5 zero 
crossings and a flight without delay has 16 zero crossings. On the contrary, for pilot P4, a flight with 
transport delay has 24 zero crossings and a flight without delay has 7 zero crossings. Using 
Lanchaster’s estimate (reference Etkin and Reid 1996, p. 172), the time period for phugoid 
oscillations at 450 kts. or 760 ft./s is 104 sec. which implies nearly 2 zero crossings in 100 seconds. It 
is not clear if the phugoid oscillations have a role here. Based on rms values of h and 8 fluctuations 
only, the performance of both pilots is nearly the same. Having fewer zero crossings in flights without 
delay as compared with the ones with delay is intuitively more acceptable. Therefore, the 
performance of pilot P4 is somewhat predictable but that of P3 is surprising. It may be interesting to 
reiterate that pilot P3 is a licensed and experienced pilot but the experience of P4 is limited to flight 
simulators only. 

Figure 4 shows the coefficient of cross correlation between aileron deflections and airplane heading 
versus time gap for the same five selected flights which have been covered in Figure 1. For four of 
the flights the graphs have similar features. For all the five flights, the largest magnitude of the 
negative correlation occurs at the time gap values between 450 ms and 725 ms. This time gap range is 
comparable with 350 to 650 ms, which is the first time gap range for 8-h correlations in Figure 1 . For 
zero time gap, an appreciable negative correlation exists in all the five flights. For both pilots, the rms 
values of aileron deflections for flights with delay are appreciably larger compared with those for 
flights without delay. The rms values of airplane heading do not have the same character as that of 
the rms values of aileron deflections. 

Figure 5 provides the coefficient of cross correlation between rudder deflection and airplane heading 
versus time gap. In the best and average flights with delay for pilot P3 appreciable changes appear in 
the correlation with time gap. In the rest of the three flights the correlation is rather insensitive to the 
time gap. The raw data in Appendix A helps to trace this insensitivity to an inactivity on the pilot’s 
part in using the rudder for heading control. The pilot P3 has used rudder movement much more 
than pilot P4. For the average flight of pilot P3 with delay, the rudder-heading correlation in Figure 5 
has similar features as the aileron-heading correlation in Figure 4. This similarity may be an indicator 
of coordinated use of rudder and ailerons. 

Figure 6 has the coefficient of cross correlation between elevator deflection and aileron deflection 
versus time gap. The correlation for the five flights does not seem to have any common features 
except that the correlation is rapidly changing between positive and negative values. 

Figure 7 has coefficient of cross correlation between elevator deflection and rudder deflection versus 
time gap. Four of the five flights show appreciable negative correlations for zero time gap, although 
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such a feature between the movements of elevator and aileron may not be worth emphasizing. In 
general, the correlation graphs for the five flights do not have any common features. 

Figure 8 presents the coefficient of cross correlation between rudder deflection and aileron deflection 
versus time gap. For different flights, at zero time gap, the different values of correlation coefficient 
appear to represent the pilot’s coordinated use of rudder and aileron. No common features are 
discemable between different flights, perhaps, because of the lack of use of the rudder by the pilots. 
The hardware facility in the present simulation system, perhaps, demotivated pilots from using the 
rudder. The system hardware does not have rudder pedals, instead, the rudder pedal movement is 
simulated by the rotation of the joystick as described in the section on “The Simulation System.” 

An overall consideration of all the correlation graphs suggests that for reasonably well defined flights, 
the cross correlations between pilot input and airplane response versus time gap reveal interesting 
features of the pilot dynamics. When the aileron-heading and the rudder-heading correlations for a 
flight are considered together, they reveal the presence or absence of coordinated use of ailerons and 
rudder by the pilot. Consideration of correlations amongst the parameters that represent longitudinal 
motion of the airplane and correlations amongst the parameters that represent lateral motion of the 
airplane support the generally used approach of separate studies of longitudinal and lateral dynamics 
of an airplane in flight. 

Additional research work may be recommended to find certain characteristic features of cross 
correlations between pilot input and airplane response. Such characteristic features would then help in 
determining different progressive stages of pilot training on simulators. 
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Table I: Comparison of the Pilot’s Performance in Their Best Flights With and Without Transport 
Delay, h’ and 8‘ are the rms values of altitude changes and elevator deflection changes, respectively, 
subscript n is used for flights with no delay and d is used for flights with 500 ms delay. 
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Figure 1: Cross Correlation Coefficient Between 5 and h 
vs. t for the Five Selected Flights of Pilots P3 and P4. B 
is the pilot’s best flight without additional transport 
delay. BD and AD are the pilot’s best and average 
flights with additional transport delay of 500 ms. (8 - 
elevator deflection, h = airplane altitude, t = time gap). 
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Figure 2: Cross Correlation Coefficient Bct ^f n 5 
LTvs. x for the Same Selected Flights Winch Am 
Covered in Figure 1. (8 = elevator deflection, 
airplane altitude). 
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figure 3: Cross Correlation Coefficient Between S and h 
vs. x for x up to 100 Seconds. For the same flights, x 
goes up to 20 seconds in Figure 1. (5 = elevator 
deflection, h = airplane altitude). 
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Figure 5: Cross Correlation Coefficient Between S, and 
y vs. t for the Same Selected Flights Which Are 
Covered in Figure 1. (8 r = rudder deflection, V = airplane 
heading). 
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Figure 6: Cross Correlation Coefficient Between 8 ^ 
5 \T x for the Same Selected Flights Which Are 
Covered in Figure 1. (8 = elevator deflection. 8. = 
aileron deflection). 
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Figure 7: Cross Correlation Coefficient Between 8 and 
5 r vs. t for the Same Selected Flights Which Are 
Covered in Figure 1. (8 = elevator deflection. 8, = redder 
deflection). 
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Figure 8: Cross Correlation Coefficient Between 8, and 
8 vs. T for the Same Selected Rights Which Are 
Covered in Figure 1. (8, = rudder deflection. 8, = aileron 
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CONCLUSION 


In the study of straight level flights on a simulator with pilot in the loop, the root mean square values 
of altitude variations for flights with a 500 ms additional delay are twice or larger compared with their 
respective values for flights with no delay. This is in agreement with the other recent studies which 
indicate deterioration in flight with 300 ms transport delay. The graphs of cross correlations between 
pilot input and airplane response versus time gap reveal interesting features of the pilot dynamics. A 
time gap of 350 to 600 ms appears to include three steps of altitude control, these are pilot’s reaction 
time, time gap between stick movement and resulting elevator deflection, and the time gap between 
elevator deflection and affected altitude change. This observation, however, requires further 
exploration because introducing a 500 ms additional transport delay does not make any visible 
impact on the correlation graphs of a flight. 

The correlation between elevator deflection and time derivative of altitude becomes zero at the same 
time gap when the correlation between elevator deflection and altitude has a local maximum or 
minimum value. An examination of the aileron-heading and the rudder-heading correlations suggest 
that similarities or lack of similarities between the two kinds of correlations correspond to the 
presence or absence of coordinated use of aileron and rudder by the pilot. Additional work is 
recommended to investigate cross correlations between pilot input and airplane response to search for 
their characteristic features which would determine different progressive stages of pilot training on 
simulators. The UNIX based simulation system that has been installed and used for the present work 
is suitable for a future study of the Pilot Induced Oscillations. The installation of rudder pedals 
suitable for the UNIX based simulator system is strongly recommended as an improvement in the 
system. 
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NOMENCLATURE 


Symbol Description 


ms 

rms 

t 

x 

H (t) 
avg H 
h(t) 
h’ 
h i + j 
¥ 

¥(t) 

A 

avg A 

8 . 

6 ’ 

8 (t) 

8, (t) 

8 r (t) 

CORR (x) 

DEP 

DAP 

DRP 

HDOT 

PSI 

SEL 

MEL 

INST 

AIAA 

NASA 

PIO 

CPU 

SGI 


Milliseconds. 

Root mean square. 

Time. 

Time gap. 

Time gap for the maximum absolute value of negative correlation. 
Instantaneous value of airplane altitude. 

Time average value of altitude for a 3 minute flight. 

H (t) - avg H, fluctuating part of altitude, 
rms value of h. 

h at the (i + j) time interval of a flight, each interval is 0.025 ms. 

Heading, zero indicates Northward direction. 

\|/ (t) - avg v|/, fluctuating part of heading. 

Elevator deflection, positive for upward deflection. 

Time average A for a 3 minute flight. 

8 at the i lh time interval of a flight, each interval = 0.025 ms. 
rms value of 8. 

A (t) - avg A, fluctuating part of elevator deflection. 

Fluctuating part of aileron deflection. 

Fluctuating part of rudder deflection. 

Coefficient of cross correlation between two signals with time gap x. 

Stick deflection for elevator or deflection of elevator, according to context. 
Stick deflection for aileron or deflection of aileron, according to context. 
Pedal movement for rudder or deflection of rudder, according to context. 
Time derivative of H. 

Heading, zero indicates Northward direction. 

Single-Engine Land (aircraft type). 

Multi-Engine Land (aircraft type). 

Instrument Rating (FAA-certification). 

American Institute of Aeronautics and Astronautics. 

National Aeronautics and Space Administration. 

Pilot Induced Oscillation. 

Central Processing Unit. 

Silicon Graphics, Inc. 


Units 

milliseconds 


instant 

seconds 

seconds 

feet 

feet 

feet 


feet 


ft. / sec. 
degree 
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APPENDIX A 


The graphs of the raw data for the five selected flights of pilots P3 and P4. These 
flights are discussed in the main text. 

B represents the pilot’s best flight without additional transport delay. 

BD represents the pilot’s best flight with additional transport delay. 

AD represents the pilot’s average flight with additional delay. 

The x-axis in all the graphs shows time from 0 to 180 seconds. The control input 
signals and the airplane response signals are plotted on the y-axis. 

H = airplane altitude (ft.). 

HDOT = dh/dt (ft./sec.). 

PSI = airplane heading (radians in this appendix, degrees in the main text). 

DEP represents joystick movement (in.) to deflect the elevator. 

DAP represents joystick movement (in.) to deflect the ailerons. 

DRP represents rudder pedal movement (in.) to deflect the rudder. 

On a proportional scale, DEP, DAP, and DRP may also be interpreted as the 
deflections of elevator, ailerons, and rudder, respectively. 


A 1 



P3 AD-dat 




Figure Al: Airplane altitude, time derivative of altitude, 
and airplane heading vs. time for the five selected flights. 
For nomenclature, please see the first page of Appendix 
A. 
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Copy of the paper “ Pilot Input and Airplane Response Cross Correlations on a Flight Simulator With 
and Without Transport Delay. ” The paper was presented by Jason Paul Williams and Yusef Ali Johnson at 
the AIAA Southeastern Regional Student Conference, held on April 10 and 11, 1997 in Atlanta, GA. The 
work reported here was performed as part of the present project under the NASA Dryden grant number 
NAG 2 - 4006. 
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Abstract 


Four volunteer pilots were required to maintain given altitude and heading on a computer based stationary flight 
simulator of the F-15E fighter airplane. One set of flights was obtained without any additional delay between the 
stick input and the appearance of airplane response. In another set of flights, a 500 ms additional delay was used. 
The deviations in altitude and heading and the fluctuations in the altitude indicated that the flights with 500 ms delay 
were appreciably deteriorated. This is in agreement with the other recent investigations. The cross-correlations 
between elevator deflection and altitude reveal interesting aspects of pilot’s performance. In particular, these 
correlations reveal that different pilots exercised a variety of techniques for maintaining straight level flights. The 
C++ source code, which is used for the computation of the desired correlations, is also included. 
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Introduction 


For human operators in a feedback control system, their ability to adapt to variations in control tasks is generally 
acknowledged as a significant advantage. Baron 1 says, “With training and proper feedback, pilots are capable of 
adapting their control strategies to a wide range of task variables .” Baron, however, reminds us that this adaptive 
capacity is not unlimited. The transport delay or a delay between input signals and their response is one of the most 
important factors that needs studies for different flight situations to understand whether they fall within the pilot’s 
adaptive capacity. The papers published by Cardullo et al 2 , Lusk et al 3 , and Middendorf et al 4 are some of the recent 
studies on the effect of simulator transport delay on pilot’s performance and on the ways of compensating delays. 
The present paper offers a study of the airplane responses to stick inputs with and without an additional transport 
delay. Four volunteer pilots conducted straight level flights on an F-15E flight simulator. The pilots used the joy 
stick to maintain a given altitude and a given heading. Two sets of flights were conducted. In one set the standard 
simulator was used and in the other a 500 ms additional delay was programmed between the stick input and the 
appearance of airplane response on the monitor. The deflections of elevator and ailerons were treated as pilot’s 
inputs, and the altitude and heading of the airplane were treated as the responses. In the present study the average and 
rms values of the elevator deflection and altitude signals as well as the cross correlations between these two signals 
are presented as the statistical tools to compare pilot’s performances. A C++ code to compute the desired cross- 
correlations from the flight data is provided in the appendix. The present paper is included in the final report on a 
flight simulation research project at Tuskegee University sponsored by NASA Dryden Flight Simulation Research 
Center under the Grant NAG2-4006, June 1994 to January 1997. 

Background 

The simulation software, which represents the dynamics of flight of an F-15E fighter, has been installed by NASA 
Dryden Flight Research Center on an SGI 4D/420 IRIS workstation. On the IRIS alone, the simulator could be 
operated by a mouse. To operate the simulator by a joystick, its graphic software has been moved to an SGI 
INDIGO 2 workstation. The two workstations are connected through ethemet. Thus the F-15E simulation software 
used in the present study operates on three central processing units, two on an SGI IRIS workstation and one on an 
SGI INDIGO 2 workstation. The pilot in the loop operation is comprised of the simulation software, the pilot, and 
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a BG System joystick called a FlyBox. Norlin 5 has provided a description of the software design, model 
development, and simulation capabilities of the system together with a description of other NASA Dryden Systems. 
The F-15E simulation is coded with a FORTRAN shell and C support routines and operates on a UNIX - based 
platform. The integration scheme in the simulation has been optimized for real time operation. The code has the 
flexibility to program a desired transport delay between the stick input and the airplane response on the monitor. In 
a simulated flight, the digital data on forty-two different flight parameters can be recorded at every 25 ms time 
interval. The flight parameters include the altitude, heading, and velocity, as well as the deflections of elevator, 
rudder, and ailerons. Provisions for turbulence and gust, although available, were not used for this study. 

Experimental Protocol 

Four volunteer pilots conducted the required flights for the experimental study on the F-15E airplane real time 
simulation. Three of the them were experienced licensed pilots and one of them had no prior experience of flying. 
In this paper, the volunteer pilots are named as PI, P2, P3, and P4. The ages, flying hours, and license status of the 
pilots are PI: 69 yrs, 3200 hrs, SEL, MEL, INST; P2: 48 yrs, 100 hrs, SEL, MEL, INST; P3: 23 yrs, 830 hrs, 
SEL, MEL, INST; and P4: 26 yrs, no flying experience. 

In every test flight a pilot was required to maintain straight level flight for three minutes duration heading toward 
North, with pre-adjusted throttle for flying at 450 knots indicated at 10,000 ft. altitude with autotrim off. During the 
flight, the pilot inputs to the joystick were limited to elevator and aileron deflections only. A measure of 
performance, named rsws score, was calculated for a flight from its altitude and heading data. An rsws score is the 
square root of the sum of the weighted squares of maximum deviations in heading and altitude. Let A h in ft. be the 
maximum absolute deviation from the required 10,000 ft. altitude. Let A 4* in degrees be the maximum absolute 
deviation from the required Northward heading. Then 

rsws = sqrt [(A^F) 2 + (A h/10) 2 ] (1) 

It is obvious that a smaller value of the rsws score represents superior pilot performance. Two kinds of test flights 
were conducted. One kind of flight simulated a normal flight without inserting additional delay between stick inputs 
and the appearance of airplane response on the computer monitor. The second kind of flight had a 500 ms transport 
delay programmed between the stick inputs and airplane response. For each of the two kinds of flights, a pilot was 



trained and evaluated by conducting several three minute flights and by calculating the respective rsws scores. An 
asymptote in performance consisted of six consecutive flights of the same kind in which the pilot did not break his 
own rsws score record. The flights of each kind that represented the asymptote in performance for every pilot formed 
the data base for further analysis. In addition to the rsws scores, other criteria were also considered to compare 
different flights. For the present paper, attention is focused on the elevator deflection, A as input signal and altitude 
of airplane, H as the response or output signal. The time average values of A and H are avg A and avg H. The 
fluctuating parts of the two signals are 8 (t) = A (t) - avg A and h (t) = H (t) - avg H, where t represents time. The 
rms value of 5 and h are denoted by 8’and h\ The cross-correlation between 8 and h for arbitrary time gap x is 
denoted as CORR (X) and defined as: 

t=T 

C ° RR(T) = TS'h' J 6(0 h(t + T) dt 

t=0 (2) 

The simulator provides data at every 25 ms interval, thus a 3 minute flight has 7201 data points for each signal. 

The cross-correlation is therefore calculated as follows: 


qwj) 

CORR(j) = — - — / 8*h 

(7200 -j) 6' h l 1+J 

1 = 1 (3) 

To cover a time gap, x, of 20 seconds for correlation, CORR (j) is calculated for j = 1,2, ..., 7201. For each kind of 
flight for every pilot, the correlation graphs, CORR (x) vs x are obtained for only three flights amongst the six 
flights that determine the asymptote in performance; these three flights are the best, the worst, and the average based 
on the rsws scores. An initial estimate of 20 second maximum time gap, x, was considered sufficient for 
correlation. Upon considerable computations we learned that longer time gaps would have been desirable. 


Fli ght Data Analysis 

In Table I, PI, P2, P3, and P4 are the four volunteer pilots who conducted two kinds of straight level flights of three 
minutes duration. The column 1 for each flight has rsws scores. A lower value of rsws indicates better pilot’s 
performance The columns 2 and 3 have rms values of elevator deflection in degrees and altitude fluctuations in ft.. 
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respectively; their lower values also indicate better pilot’s performance. From Table I it is apparent that 500 ms 
transport delay has caused deteriorated pilot’s performance. Middendorf et al 4 has noted that 300 ms delay resulted in 
significant deterioration of pilot’s performance. Lusk et al 3 have conducted an elaborate study on the delay 
compensation. With 300 ms delay they observed more significant deterioration in heading performance than in the 
altitude performance. Consideration of pilot P2 in Table I brings up an interesting aspect of pilot’s strategy. Unlike 
other pilot’s, the pilot P2 has smaller rms value of elevator deflection for the flight with transport delay than the 
flight with no transport delay. To obtain improved understanding of pilot’s performance, cross-correlations are 
obtained between the elevator deflection and altitude. From amongst the flights that represented asymptotic 
performance of a pilot, only the best, the worst, and the average, are processed to report the cross-correlations. Thus 
for the two kinds of flights by every pilot, a total of twenty-four cross-correlation curves are reported. Figures 1, 2, 
3, and 4 show the cross-correlations between the fluctuating parts of elevator deflection, 5(t) and the altitude h(t + x) 
for the time gap, x ranging between 0 and 20 seconds. The positive value of 8 represents upward deflection of the 
elevator that gives decrease in lift on the tail and which should result in nose up tendency of the airplane. Thus the 
positive 5 is expected to result in positive h. Twenty-three out of twenty-four correlation graphs show negative 
correlations for the time gap values near zero. While conducting their test flights, the pilots had observed that in the 
event of leaving the stick free, the airplane had a tendency to climb up. To arrest the climb up tendency, the pilot 
moves the elevator down that gives rise to negative correlation. Thus the negative values in the first part of the 
correlation curves reveal lack of trim. Let x m be the time gap at which the negative correlation has the maximum 
absolute value. The x m is ranging between zero and 600 ms and it does not seem affected by the transport delay. We 
would expect x m to provide a measure of the pilot’ s reaction time delay but its apparently unpredictable variations do 
not meet the expectations. It may be noted here that Roskam 6 assumes 100 ms reaction time for test pilots and 120 
to 200 ms for other pilots. For larger time gaps, most correlation curves show appreciable positive values. It is 
suspected that they are caused by the Phugoid oscillation of the airplane, reinforcing the desire to extend time gap 
delays to larger values. 
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Conclusion 


For all pilots observed in the study during straight level flight on a simulator, the additional transport delay of 500 
ms, between the stick input and appearance of airplane response on the monitor, results in appreciable deterioration 
of the flight. This is in agreement with the other recent studies which indicate deterioration in flight with 300 ms 
transport delay. The cross-correlations between elevator deflection and altitude signals support a belief that in a non- 
perturbed flight, different pilots exercise a variety of techniques for maintaining straight level flight. It is proposed 
that cross-correlations between the elevator deflections and aileron deflections shall provide an improved 
understanding of pilot’s technique on a flight. 
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Nomenclature 


Symbol 

t 

T 

A 

avg A 

8 

5 ’ 

5 (t) 

h i+j 

h’ 

h(t) 

H 
AH 
avg H 
V 

A *¥ 
rsws 

CORR(x) 

SEL 

MEL 

INST 


Description 
Time 
Time gap 

Time gap for the maximum absolute value of negative correlation 
Elevator deflection, positive for upward deflection. 

Time average of heading for a 3 minute flight. 

8 at the time interval of a flight, each interval is 0.025 ms. 
rms value of 8. 

A (t) - avg A, fluctuating part of heading. 

at (i + j) the time interval of a flight, each interval is 0.025 ms. 

rms value of h. 

H (t) - avg H, fluctuating part of altitude. 

Right altitude. 

Maximum absolute deviation in altitude, from 10,000 ft. 

Time average value of altitude for a 3 minute flight, ft. 

Heading, zero for Northward heading. 

Maximum absolute deviation in heading, from zero deg. 

Square root of the sum of the weighted squares of maximum 
deviations in heading and altitude. 

Coefficient of cross-correlation 8(t) and h (t + x). 

Single-engine Land (aircraft type) 

Multi-engine land (aircraft type) 

Instrument Rating (FAA certification) 


Units 

instant 

seconds 

seconds 

degree 

degree 

milliseconds 

degree 

milliseconds 


feet 

feet 

feet 

feet 

degree 

degree 
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Pilot 

Flight Without Transport Delay 
1 2 3 

Flight With Transport Delay of 500 ms 

1 2 3 

PI 

26.4 

0.60 

73.8 

52.2 

1.02 

77.1 

P2 

29.3 

0.51 

51.2 

33.5 

0.20 

66.4 

P3 

3.2 

0.23 

mm 

5.8 

0.45 

18.2 

P4 

3.7 

0.25 

mm 

6.3 

0.72 

18.7 


Table I: Comparison of Pilot’s Performance in Flights With and Without Transport Delay. 

1 - Average rsws for the best 5 flights, 2 - 8‘ deg in the best flight, 3 - h’ ft in the best flight. 
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Pl_B.bdep 


*P1 B.bdep ' 

0 4 :OKtia« .0.193746 • tau -15.35 
COBft-la --0. 366041 • cau -0.45 
iv«l - 9947.47 av«OI? - >2.05784 
rail « 73.8086 vmm0€9 - 0.603538 
0.2 ^ 


Pt_BD.hdep 


' PI BD.hdep 

0 B 8 »ax -0.113324 • tan -7.4 - ‘ 

OKBaia ->0.509733 8 tau -0.3 
v*B - 9969.2 tvaOEP - >1.62183 
- 77.0513 raaOKP • 1.01876 


10 15 

TAU 


Pl_A.hdep 


O 4 IpOKftaaH -0.294876 • tan -11.55 
U ‘ COBBaln --0. 372214 • tan -0.35 


1 Pt_A.hdep 


10026.3 «v«0£2 

126.69 ryW 


a^_ >1^7 9594 


Pl_AD.hdep 

# [ ’ * P^AD.hdep' 

:otta«i -0.363204 6 tu -19.975 
6 COIlal* --0. 62468 • can -0 

avail - 9603.67 avsOSP - >0.692137 
.4 rail - 1035.78 rasOSP - 1.32935 


Pl_W.hdep 

°-*r " Pl_W.hdep ’ 

0.3 COttaaa -0.266918 • t.u -19.975 
COBBala ->0. 0736797 • tau - 
aval ■ 10123.6 avaOIP ■^iii llli I l 
°* 2 ras6 - 204.372 r sjJfcS^-Co . 7 5 9 1 4 


PI WD.hdep 


•PI WD.hdep 

Ollaat -0.263664 t tan -19.973 
Qllaia ->0.545517 • tau -0 
vs« « 9953.01 avaOIP - 0.220115 
..A - 16S0.04 raaOIP • 2.29504 


Figure 1: Cross -Correlations Between Elevator Deflection and Altitude of the Airplane for Pilot PI 

B. A, and W are the pilot's best, average, and worst flights without additional transport delay 

BD, AD, and WD are the pilot's best, average, and worn flights with additional transport delay of 500 ms. 
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Figure 2: Cross-Correlations Between Elevator Deflection and Altitude of the Airplane for Pilot P2 

B. A. and W are the pilot’s best, average, and worst flights without additional transport delay 

BD. AD. and WD are the pilot’s best, average, and worn fliah* with additional transport delay of 500 ms. 
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Figure 3: Cross -Correlations Between Elevator Deflection and Altitude of the Airplane for Pilot P3 

B. A, and W axe the pilot's best, average, and worst flights without additional transport delay 

BD, AD, and WO are the pilot's best, average, and worst flights with additional transport delay of 500 ms. 
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Figure 4: Cross-Correlations Between Elevator Deflection and Altitude of the Airplane for Pilot P4 

B. A, and W are the pilot’s best, average, and worst flights without additional transport delay 

BD. AD. and WD are the pilot's best, average, and worst flights with additio n al transport delay of 500 ms. 
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Appendix 


The complete program written in C++ is enclosed. The program reads the F-15E simulation data after the data file 
has been converted from binary to ASCII. The program was compiled on IRIX 5.3 C/C++ compiler 4.0.2 on the 
SGI Indigo2 workstation. Pages 15 to 17 contain the Correlation Source Code. Pages 18 and 19 have the 
Correlation Header. Page 20 has the Script Source Code. The program prints out the calculated correlation values 
and the jargon for plotting the graphs on a GNU Plot Program. The Script file separates the jargon from the data, 
and it creates a script file for the GNU Plot Program. 
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♦ include < f stream. h> 

♦ Include <string.h> 

♦ include <math.fa> 



♦ include 0 corr 1 . hpp * 

/* 

correlation header */ 

— 

♦define 1 7200 

void main ( ) 

( 

char fileName (80 ] ; 

/* 

maxi mu, number or points */ 


/* 

input filename */ 



char FileName( 30 ] ; 

/* 

output filename */ 


char buffer C 801 ; 

/* 

buffer */ 


char caxt(255] ; 

/* 

notes character variable */ 

— 

char CAtlOI ; 
char CB ( 10 j ; 
char FileOataCSOT ; 
char FileloadC801 ; 

/• 

name of A(t ♦ tau) */ 

/• name of B(t) V 


int chan; 

int i, ii > 0, j * 0,k » 0; 

/• 

number of channels */ 

— 

int nC101 ; 
int w; 
int cf; 

/•data pointer */ 

/* correlation type •/ 

/• maxium correlation factor 


int point; 
double U; 

/* 

enter maxima number of point */ 



double d; 

/• 

“tau value at maxium correlation " 


double e; 

/* 

tau value at minium correlation * 


double Numb; 

/* 

data buffer */ 


double yu; 

/• 

maxium correlation */ 


double yl; 

/• 

minium correlation */ 

— 

double xznin; 

/• 

minium x value of graph •/ 


double ymin; 

/* 

minium y value of graph */ 


double xmax; 

/* 

maxima x value of graph */ 

— 

double ymax; 

CORR DATA, Scorr, corr; 

/* 

maxima y value of graph •/ 


rename: 

eput << * Filename: •; 

cin » filename; /* filename of the ascii */ 

cout « * Enter the maxium correlation factor * ; 

cin >> cf; /* enter maxium correlation factor * f 

cout << • Eater number of points *; 

cin >> point; /* enter m**™ number of point */ 

cout << * Please choose and correlation (HJDE?*1» PSI_0AP*2, PSI_ DRP«3) : •? 
cin » wr /• correlation selection •/ 

/* opens the data file, if not file exist then the file name is asked 

for again */ 

if stream f ia( filaName, ios : :nocreace) ; 
if (!fin) 

( 

cout « • Unable to open :: • << filename << * *** file not found *** \n\n 0 ; 

cout << 0 Please re-enter the filename. \n\n 0 ? 
goto rename; 

) 

cout << 0 \n Filename 0 << fileName « 0 was found. \n* ; 

/* reads unneeded j&rgen left by the getOata p r ogr a m */ 
for ( i * 0; i <5; i+e ) 

( 

fin » buffer; 

> 

fin » chan; /• reads the number of Chanels */ 

cout « *\n This file has 0 << chan « 0 channels. \n*; 

/* reads and sorts the number of channels */ 
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for ( i* 0; i < (rhan»2) ; i++ ) 

( 

fin >> buffer; 

if ( ! (stranp (buff er, *T* ) ) ) ( n[Q] * j;) 

if ( ! (strcmp (buf f or, *H* } } ) ( n[l] * j?} 
if ( » (strcmp(buf for. *PSI* J ) ) ( n[2] •» j?) 
if (! (acrcap (buffer, * DAP* )) ) ( n(3] * j?} 
if (» (strcmptbuffer. *DEP*) )) ( n(41 * j;> 
if ( 1 (strcmp (buffer, *DRP*) ) ) { n(51 ■ j?) 

j*** 

} 


/* reeds and sores the data to thier correct arrays */ 
while ((fin >> Numb) && Oc < (point)) ] 

( DATA.TCk] = NUmb; ) 

( DATA.Ktkl * Numb;) 

{ DATA.PSKk) * Numb? > 

( DATA. DAP Ckl = NUmb? > 

{ DATA. DEP (k) * Numb;) 
(DATA.DRP(k) * Numb? ) 

{ii * 0; k*+? ) 


/* closes the input file */ 

to the corrolation header*/ 

DATA.cfscf ? 

DAT A. ww i 
DATA.k-k? 

DATA.corr ( ) ; 


if ( ii ** n(0] ) 

if ( ii *■ n(l] ) 

if ( ii ** n(2] } 

if ( ii ** n(3J ) 

if ( ii ** n(4] ) 

if ( ii *- a(5] ) 

Ll++t 

if (ii*- (chan+I) ) 

> 

fin. close ( ) ? 

/ * sets values initalize 


/* finds the «*w4t«n upper and lower points of the correlation / 
yl * DATA-CORPAt G 1 ;d*0;yu » DATA.CORPA(0 1 ; e-0 ; 

for (i«0?i<cf? i*+ ) 

* if (yl > ( DATA . CORPA ( i J ) ) (yl * (DATA. CORPAdl ) ;d-i # ( .025) ? > 

> 

for ( i * 1; i < cf? i** ) 

( if(yu < (DATA.CORPACil) ) (yu » { DATA. COSPAC ill ;«-i* ( 025) > 

if(fabsiyu) < fabe (yi) ) {U*yl? } cout <<tJ? 

if (fabs(yu) > fabs (yl) ) (tJ*yu? ) cout «U; 

xzsax * ( (DATA.TAUCcf-H ) # .2fDATA.TAa(cf-ll ) ; 

vwi n m 0 ? 

ymin * — { f abs (C>D* . 50 ) ) ; 
yraax * ( fabs ( U>U*.50J); 


Rename: 

cout << * Filename: *; 
cin >> FileName; 

cout << * Enter one line of co nanen ts 
cin. ignore ( 1, * \n * ) ? 
cin.getline ( text, 253) ? 


/* 

for data file 
/• 


name of the output file */ 

:\n*? .. u ./ 

enter notes on the graph / 


/* 


sorts the the file name */ 

• strepy (FileBata, FileName ) ; 
strepy (Fileload, FileName) ? 
street (FileLoad, * .gif) ; 
if (w~l) 

street (FllaData, * .hdap* ) ; 
street (CA, *H*>? 
street (CB, *DEP* ) ? 
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} 

if (w»*2 ) 

( 

strcatt (FilaData, • .pdap - ) ; 
s treat (CA. • PSX* ) ; 
streat (CB. • DAP - ); 

> 

if (w»«3 ) 

( 

streat (FilaData, * .pdrp* ) ; 
streat (CA, • PSI* ) ; 
streat (CB.*CRP*J; 

) 

if (w»4 ) 

( 

streat (FilaData. *. edrp* ) ; 
streat (CA* *DSF* ) ; 
streat ( CB , * DRP ■ ) ; 

) 

if (w**5) 

{ 

streat (FilaData. ■ . edap * ) ; 
streat (CA, ■DEF* ) ; 
streat (CB.-DAP*); 

> 

/* opans output fila •/ 

of stream f out (FilaData* ios : : no rep lac a ) ; 
if (If out) 

cout « * Unable to craata and write to * « Filename « 
cout << * Plaasa re— enter the filename. \a\n* ; 
goto Rename ? 

> 

/* writes acripe jargen */ 

foutt « *t set title \ ' • «FileOata« * \ ' \n* j 
fout « •# set xlabel \ *TAU\ ' *«* \a* ? 
foutt « see ylabel VCORRV 0,-10*«'\n-; 

fout « •• see label 1 VCORRmax «*«yu«* % tau **«e«*\ at « 

foutt « •* see label 2 \ 'CORRmin -*«yl«’ 3 tau e*«d«*\' at << 

foutt « •# see label 3 Vave*«CA«* * • <<OATA. aveA« * av * #<< “ <K 

fout « •# set label 4 \'rms*«CA«* * • «DATA. rmsA« * nas*«CB« 

fout << *• plot (x»*<< xmin <<•:•<< xm a x ; 

fout <<•] [•« ytnin <<• : *«y»ax«* 1 V •«FileData«*\ # \a* ; 

fout <<•# tau(i] corr(i] “ « *\a*; 

fout <<*# • << fileName << - \a*j 

fout <<•* * << text << *\n*; 
fout <<*• tau(il corrHCi] " << 

/* writes results to output file */ 
for ( i » 0? i < cf; i*+ ) 

fout << DATA . TAU ( 1 ] «*\t* << DATA . CORP A ( i 1 « 9 \ti m ; 

) 

f out. close O ? 

cout « •\nFilename •« FilaData « * was created. \n ; 
cout << * \nNow run script and use ■<<FilaData«* script 


/* closes 

file.\n* ; 


file alreadv 


xmin*. 9 «*, * 
xmin*. 9 «•, * 
=* # «DATA.ax 
1 = *«0ATA.rs 


output file * 
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♦ include <math.h> 
•include <fstream.h> 


•define 1 7200 

class CORK 

{ 

public : 


int i, j, Jc? 
int w; 
int cf; 

toubl.Ttll. Htll. PSItU. DAP (11. DEP ( 1 ] > DRPtll. TA0(U . CORPAU1 
double t, h, psi, dap, dep, drp; 


double aveA, avefl ; 
double sqdA, sqdB; 
double nsnsqA, mnsqfl ; 
double rasA, rmsB; 
double corp# corn; 
double saA# sbB; 
double sasqA# sasqB ; 
double smcorn, amcorp; 
double A(l] ; 
double BC1] ; 
double di£S[l]; 
double difAtl] ; 
double corpA; 
double cornA; 
double tau; 


/ 


coniine > 

-CORR( ) ( } 
void Scorr () 

* T(z]*t, HC*]*h. PSI(zl 3 psi. DAP ( x ] *dap • DEPtzl*d«P» DRP(zl*drp; 


) 

void cost 0 


secs all sum variables to zero 

places correct signals to be corelated to 

sums A and B */ 

■ saA * 0, smB » 0, sasqA ■ 0, sasqB * 0, 


A and B 

* 0, 


sacorp 


0 


if (w**L) 

for ( i * 0? i < k; i+«* ) 

{ Ati] * H[i] ; B C 1 1 3 DEPCil; 
satA 3 smA ♦ A(i] ; 
smB * smB ♦ B[i] ; 


) 

) 

if (w»*2) 

for ( i * 0; i < k; i+* ) 

( ACi] * PSIti] ?B(i] » DRPtil; 
saA 3 saA ♦ A(i] ; 
smB * saB ♦- 8(11 ; 


) 

> 

if (w»3) 

{ 
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DAPti] ; 


for ( i « 0 ; i < k; i++ } 

( 

A(i] « PSI(i];BCi] » 

3 smA «■ A(i] ; 
im B ♦ B(i] ; 

) 

> 

if (w**4) 

{ 

for ( i * 0; i < k; ) 

ACil * DRPCil ;B[i] * DE?(i] ; 
smA * smA +■ A[i] ; 
saB m sad * B(il; 

) 

} 

if (W”5 ) 

( 

for ( i * 0; i < k; i*+ ) 

1 Ati] a DAP(i];B(i] * DEPCi]; 
a smA + A(i] ; 
i smS ♦ B [ i 1 ; 

) 

> 

aveA * smA / (k) ; ' 

&v«B « saB/ (JO ; /* 

for (i * 0; i < k;i+«*) 

( dif B [ i ] » <s Cil - aveB) ; /* 

difACi] * (A(il - aveA); /* 

sqdB * pow<difB(i] # 2.0); 

•qdA a pow< difACi] . 2.0); 
smsqA a smsqA «■ sqdA; 
mqB ■ smsqB ♦ sqdB; 

> 

mnsqA * sasqA/ Oc) ; 
masqB * smsqB/ Oc) ; 

rasB a (sqrt (mnsqfl) ) ; /* 

rmsA * (sqrt (mnsqA) } ; / f 

for (j ■ 0; j < cf ? j**) 

* for (i a 0; i < (k-j+1) ; i+«-> 

( corp a cLifBCi) ♦ difACi* j 1 ; 
saeorp * sncorp + corp; 

} 

C&u « j * 0.025; 

corpA a saeorp/ ( (k - j) * rasB * 
saeorp a 0 ; 

TAUCjl * tau; 

CORPA C j 1 = corpA; / 

) 


avarags of A */ 
average of B */ 


dif f srencs between BCU “d aveB */ 
diffaranca batwaan B[il and avaB / 


ras valua of A * / 
ras valua of B */ 


ms A) ; 


/* correlation coeficient 


correlation coeficient array */ 


•/ 
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buffer[255] , WRCTEUSS]** • , f ileLoad[80I . scrl [30 ] 


* was found. \n" 


• include <f stream. fa> 

•include < string. ii> 
void main ( } 

( 

char fileName (80] , 
int i.n; 

// enter filename 

cout << * Filename: *; 

cin » fileName; 

cout « ■ \n Filename * « fileName << 
ifstream f in ( fileName. ios : :aocreare> j 
strlCl] » '\0'? 
n ■ 0; 

while < fileName Cnl 

( 

strl(0] » fileName (n] ; 
street ( fileLoad, strl) ;n++? 

) 

s treat ( fileLoad, • .gif *) ; x 

cout << * \n Filename * « fileLoad << * was created. \n\n 

cout << " Now run gnuplot and load *« fileLoad<< \ ' 

of a cream t ouc ( f ilaLoad. ios : :nor«plac«) t 
for 

fin.g«e(buff«r. 2) ; fia.goclino(WRITZ.255) ; 

fouc << WRIT£<<* \n" ; 


) 
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